function y=pbj(size, power)
    noise = awgn(zeros(1, size), power);
    
    %num_samples=16000;
    %noise=(randn(num_samples,1)+j*randn(num_samples,1))/sqrt(2);
    bw=1/20; % bw=1 is sampling rate = 20 MHz
    [b,a]=cheby2(11,70,bw); 
    filtered_noise=filter(b,a,noise);
    fc=5.5/20;  % fc=1 is equivalent to sampling rate = 20 MHz
    carrier=exp(1j*2*pi*fc*(0:size-1));
    y=filtered_noise.*carrier;
    
%     center = 5.5e6;
%     F = [4.8e6 5e6 6e6 6.2e6];
%     A = [0 1 0];
%     DEV = [0.001 1 0.001];
%     [N,Wn,BTA,FILTYPE] = kaiserord(F,A,DEV,20e6);
%     b = fir1(N, Wn, FILTYPE, kaiser( N+1,BTA ),'noscale');
%     a = [1];
%     
%     y = filter(b,a,noise);
end